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We consider hexagonal cellular automata with immediate cell neighbourhood and three 
cell-states. Every cell calculates its next state depending on the integral representation 
of states in its neighbourhood, i.e. how many neighbours are in each one state. We 
employ evolutionary algorithms to breed local transition functions that support mobile 
localizations (gliders), and characterize sets of the functions selected in terms of quasi- 
chemical systems. Analysis of the set of functions evolved allows to speculate that mobile 
localizations are likely to emerge in the quasi-chemical systems with limited diffusion 
of one reagent, a small number of molecules is required for amplification of travelling 
localizations, and reactions leading to stationary localizations involve relatively equal 
amount of quasi-chemical species. Techniques developed can be applied in cascading 
signals in nature-inspired spatially extended computing devices, and phenomenological 
studies and classification of non-linear discrete systems. 
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1. Introduction 

Localisations — compact and long-living local disturbances of a medium's charac- 
teristics — are becoming hot topic of interdisciplinary non-linear sciences 
They can be found in almost any type of spatially extended non-linear systems, 
from liquid crystals to monomolecular arrays to reaction-diffusion chemical me- 
dia ffl. From a computer science point of view, localizations are ideal candidates for 
elementary processing units in 'free-space0, or collision-based^ computing devices. 
The collision-based computation is rooted in logical universality of Conway's Game 
of Life Fredkin-Toffoli's conservative logic and Margolus's physics of com- 
putation 1121, In collision-based computing, quanta of information are represented 



^Tcrm coined by Jonathan Mills, USA 
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by compact patterns traveling in an 'empty' space and performing computation 
by mutual collisions. The absence or presence, and type, of traveling patterns en- 
code values of logical variables. The trajectories of patterns approaching a collision 
site represent input variables. Trajectories of the patterns ejected from a collision 
represent the results of logical operations, i.e. output variables. 

Amongst many possible natural systems supporting localizations, reaction- 
diffusion chemical media seem to be most appropriate candidates for experimental 
implementation of collision-based computing architectures. The medium does not 
require expensive or over-sophisticated equipment, can be handled and tuned with 
relative ease, and allows for visual identification of experimental results. For exam- 
ple, there is a particular type of reaction-dif fusi on chemical system, the Belousov- 
Zhabotinsky reaction in sub-excitable mode that supports the existence of lo- 
calized wave- fragments (somewhat analogous to dissipative solitons which can 
play the role of the 'billiard-balls' in a collision-based computing system EEl 

In our previous works we have designed and studied a range of hexagonal 
cellular-automaton models of reaction-diffusion excitable chemical systems, par- 
ticularly those with concentration dependent inhibition of the activator I 4 | 21 | 5 | 
have analyzed three-state totalistic cellular automata on a two-dimensional lattice 
with hexagonal tiling, and discovered a set of specific rules that support a variety 
of mobile (gliders) and stationary (eaters) localizations, and generators of local- 
izations (glider guns). In 1^ we demonstrated that rich spatio-temporal dynamics 
of interacting localizations and generators of localizations can be used in imple- 
menting purposeful computation, including signal routing, multiple-valued logical 
operations and finite state ma chines. Despite success of preliminary studies, and 
some techniques developed to pinpoint 'best' rules supporting localizations, we 
remained somewhat puzzled and uncertain on whether rules manually selected are 
good representatives of a set of localization-supporting cell-state transition rules. 
We therefore applied the full power of evolutionary computation methods to evolve 
and select all possible rules that support mobile localizations in two-dimensional 
hexagonal ternary state cellular automata. Results of these studies are discussed in 
present paper. 

The paper is structured as follows. In Sect. [2] we introduce and define ternary 
state hexagonal cellular automaton. The evolutionary algorithm used to breed 
glider-supporting rules is outlined in Sect. [3] The evolved set of glider-supporting 
rules is characterised in Sect. HI In Sect. Owe present a set of quasi-chemical reac- 
tions derived from the evolved rules. We demonstrate, in Sect.[6l that by integrating 
the most common glider-supporting rules in one set, we select the rules support- 
ing only stationary localizations. Outcomes of the present studies are discussed in 
Sect. [71 
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2. Reaction-diffusion hexagonal cellular automaton 

We study a totalistic cellular automaton, where a cell updates its state depending 
on just the numbers, not positions, of different cell-states in its neighborhoods. We 
consider a ternary state automaton. One cell-state, S, is a dedicated substrate state: 
a cell in state S, whose neighbourhood is filled only with states S, does not change 
its state {S is a an analogue of quiescent state in cellular automaton models). Two 
other states, A and B, are assigned to be reactants. 

The cell-state transition rule can be written as follows: 

x'+' = f{aA{xy,aB{x)\as{xy), 

where (Jp{xY is the number of cell x's neighbors with cell-state p E {A,B,S} at 
time step t. As for all classical cellular automata, all cells updates their states 
synchronously in discrete time-steps. Our automata are based on a two-dimensional 
lattice with hexagonal tiling. The neighborhood size is seven: the central cell and 
its six closest neighbors. 

To give a compact representation of the cell-state transition rule, we adopt the 
formalism inSl^ and represent the cell-state transition rule as a matrix M = (Mij), 
where < z < j < 7, < i -I- j < 7, and My e {A, B, S}. The output state of each 
neighborhood is given by the row-index i (the number of neighbors in cell-state A) 
and column- index j (the number of neighbors in cell-state B). We do not have to 
count the number of neighbors in cell-state S, because it is given by 7 — (i -I- j). A 
cell with a neighborhood represented by indexes i and j will update to cell-state 
Alij which can be read off the matrix. In terms of the cell-state transition function 
this can be presented as follows: x*^^ = M^^(^^y„^^r^y . 

How do matrix M entries correspond to phenomena in reaction-diffusion chem- 
ical systems? The entries Mio = A and AIqj ^ B, i > 0, symbolize the diffusion 
of reactants A and B. Entries Alij = p e {A,B,S}, can be interpreted as a 
quasi-chemical reaction 

iA + jB — > p, 

where i molecules of species A react with j molecules of sp ecies B to produce species 
p. See detailed interpretations of some particular cases in I '^ l ^^ l ^ l. 

3. Breeding glider-supporting rules 

W e have emp loyed evolutionary computation techniques developed by Sapin et 
al I 13 | 14 | 15 116| fgj. evolving cellular automata which support mobile localizations 
(gliders). We used an evolutionary algorithm that incorporates aspects of natural 
selection or survival of the fittest. It maintains a population of structures (usually 
initially generated at random) that evolves according to rules of selection, recom- 
bination, mutation, and survival, referred to as genetic operators. A shared 'envi- 
ronment' is used to determine the fitness or performance of each individual in the 
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population. The fittest individuals are more likely to be selected for reproduction 
through recombination and mutation. 

The search space is the set of cell-state transition rules (each rule represents 
a unique cellular automaton). An automaton of this space can be described by 
determining what will become of a cell in the next generation, depending on its 
neighbours. An individual is an automaton coded as a bit string of representing 
the values of a cell at the next generation for each neighbourhood state. A string 
is composed of two sub-strings. The first substring represents the neighbourhood 
states used by mobile localizations (gliders) and their values are determined by the 
evolution of the glider. The second substring is initialized at random. The search 
space contains 3^^ possible cell-state transition rules. 

A fitness function was computed as follows. Random configurations of cells are 
evolved by the tested automaton. After this evolution, the presence of gliders is 
checked by scanning the result of the configuration of the cells. The value of the 
fitness function is the number of gliders that appeared divided by the total number of 
cells. The 36 bits of each individual are initialized at random. The mutation function 
consists of mutating one bit among 36, while the recombination is a random point 
crossover. An elitist strategy in which the best individual of population is kept 
is used. The value of the fitness function and the generation of the best rule are 
memorized. If after ten new generations the algorithm has not found a better rule 
the algorithm stops. 

4. Likehood of gliders 

From the set of evolved functions we calculatecQ a set of matrices = (^y)i 
where i is the number of neighbours in state A, and j is the number of neighbours 
in state B, z G {A, B, S}, which reflect likehoodness for any particular transition to 
contribute to glider dynamics. If one arbitrarily selects a glider from the set of thou- 
sand gliders discovered, then the local transition function supporting the glider have 
transition z amongst its state of transitions with probability F^. We have 

also computed matrix F"^ which indicate likehood of any particular state transition 
to be redundant, i.e. not necessary for supporting gliders. Exact structure of the 
matrices is shown in Fig. [TJ The matrices provide an adequate guiding choice when 
experimenting with randomly selected local transition functions, see demonstration 
applet inl^. 

Just a few examples of mobile localizations discovered are shown in Fig. [21 Most 
typical localizations are gliders, which change their configuration every second time 
step (Fig. [2^b), and puffer trains — gliders leaving trail of breathers behind them 
(Fig. [2];d) . To complete this section we would like to stress that the evolutionary al- 
gorithm employed was tuned to select transition functions supporting localizations, 

''The probabilistic matrices are calculated by Emmanuel Sapin, in the framework of EPSRC project 
EP/E005241/1 
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Fig. 1. Glidcr-likohood matrices. 
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Fig. 2. Examples of localizations observed in rules presented by glider-likehood matrices. State 
A is shown by circle, state B by solid disc, (a) glider propagating East; (b) glider propagating 
North-East; (c) puffer train propagating North-East; (d) puffer train propagating East. 



not selecting rules with manageable dynamics of localizations. Therefore, those who 
will start experimenting with matrices presented in Fig. [1] may be disappointed to 
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Fig. 3. An example of a typical configuration generated by glider-supporting rules. Gliders mul- 
tiply during collisions, so the whole space is filled by visually disorganized but in reality highly 
ordered pattern. 

see that, in most cases, development of an automaton from initial random configu- 
rations leads to disorderly looking configurations (even if the patch of initial stimu- 
lation was small enough). This is because gliders inhabit such spaces in abundance, 
they interact one with another, produce more gliders in result of their interaction, 
and populations of swarming gliders look like quasi-chaotic patterns for naked eyes 

(Fig. ED. 

5. Quasi-chemical reaction 

In order to characterize sets of possible quasi-chemical reactions, which lead to 
the formation of stationary and mobile perturbations of reagent concentrations, we 
provide isolinc representation of the glider likehood matrices in Fig.[4l Based on the 
distributions shown in Fig. |4] we can speculate that 

Finding 1. A probability, or rate, of dissociation of reactants decreases with the 
increase of number of the reactant molecules in each one local site of the medium. 

See topology of isolines in Fig. [4^. 

Finding 2. Reactions leading to production of reactants A and B are most likely 
to involve not more then four (usually between one and three) molecules of the 
reactants. 

Distributions shown in Fig. [4^ and b reach their ground-zero levels (no reaction 
take place) for the sites where satisfying the condition: Sum of molecules A and 
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molecules B does not exceed four. 



Finding 3. Reactant B is more likely to be produced during quasi-chemical reactions 
derived from F-matriccs. 



See Fig. [4^ and b to compare area of elevations. 
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Fig. 5. Concentration dynamics of A, B and S species in a well-stirred reactor. Simulated on CKS 
(Chemical Kinetics Simulator, IBM) with the following parameters: pressure and temperature are 
kept constant, volume is not tracked. lOOK particles were simulated. Initial concentration of each 
species is 1 mole/1. 

Basing on these findings we generalize glider likehood matrices to a set of ab- 
stract, quasi-chemical, reactions as follows: 
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A typical scenario of how the system (1) behaves in a well-stirred reactor is shown 
in Fig. [51 We have confirmed in the computational experiments that the reaction 
scheme developed represents an oscillatory chemical system, where concentration 
of substrate is significantly higher than concentrations of reactants A and B. This 
indeed conforms with the nature of spreading localizations and pulsating behavior 
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of generators of loealizations, or glider guns. The presence of both types of reactants 
in a locahzation is essential for the majority of propagating localisations, which is 
identified by the similarity of the dynamic of concentration profiles of reactants A 
and B shown in Fig. [5l 



6. Reductions of transitions functions 

What would be the behaviour of automata governed by the most likely cell-state 
transition rules? In order to find an answer we decided to reduce matrices of glider 
likehoods to their most articulated forms. First, we have made matrices more sym- 
metrical by the following procedure. If — < e then we updated entries {ij) 
and [ji) as follows: F-^ = F-'^ = [^{Flj + Fj'^)\ . We do this because we assume that 
in evolutionary experiments reactants A and B have the same diffusion coefficients. 

Then we perform the following operation. If \F.jj — F^\ > 0.2 and \F^^ — Ffj \ > 
0.2, or i^^ = and Ffj = then = 1. We also discard entries of (0,5) and 
(1,4) of the matrix as being negligibly small, and omit the matrix F"^ out of 
consideration. After these operation, we got the following matrix 

" [0,1,2] 0' 
1 [0,1,2] 
[0,1,2] [0,1] 
^ [0,1,2] [0,1] 
[0,1] 




where Rij is a set of all possible states that a cell x can take if its neighbourhood 
u{x) has i cells in state A and j cells in state B. We foimd that 



Finding 4. Starting its development from a random configuration a cellular automa- 
ton governed by a local transition functions expressed by matrix R exhibits still and 
breathing localizations. 

Most common localizations observed in cellular automata with local transition 
rules from the class R are shown in Fig. [6l Minimal weight localizations are still 
triples of cells in state A or B, or triples switching between A and B (Fig.[6K). These 
are followed by compact still patterns of four cells in A or B states, see examples 
in Fig. [BJdc, switching compact patterns of four non-resting states Fig. |6li; and, 
wish-bone shaped clusters of A and B states Fig. |6j2-h. The next heavier common 
localizations are the switching patterns of six Fig. [6j and eighth Fig. [6j non-resting 
states, and, the breathing localization Fig. [6l< oscillating between 12 and 13 non- 
resting states. 
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Fig. 6. Stationary localizations typically observed in automata governed by cell-state transition 
rules from class R. State A is shown by a circle, state B by a solid disc. 



7. Conclusion 

We have hybridized the paradigm of reaction-diffusion cellul a.r automata l ^ l ^-^ l ^ lwith 
evolutionary techniques of breeding glider-supporting rules I 13 | 14 | 15 116j statisti- 
cally evaluate the set of all possible totalistic cell-state transition functions which 
support mobile and stationary localizations. We calculated exact structures of glider 
likchood matrices and interpreted them in terms of abstract reactions. We demon- 
strated that quasi-chemical systems derived from glider-supporting rules exhibit 
classical dynamics of excitable chemical systems. We obtained computational ex- 
periment evidences that by reducing glider likchood matrices to their strong com- 
ponents we obtain a set of local transition rules that exhibit exclusively stationary 
localizations. Results of the research undertaken provide a priceless tool for design- 
ing collision-based computing schemes in spatially extended non-linear systems. 
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